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ABSTRACT 

This  report  contains  sixty  maps  presenting  computational  data 
pertaining  to  the  geometric  properties  of  the  geomagnetic  field 
and  a  brief  related  theoretical  description  to  make  this  material 
more  meaningful.  The  maps  display  (1)  the  values  of  total  field 
intensity  (B)  and  Mcllwain's  magnetic  shell  parameter  (L)  for 
altitudes  from  zero  to  two  thousand  kilometers,  (2)  the  detailed 
value  of  total  field  intensity  near  the  South  Atlantic  anomaly, 

(3)  the  traces  in  both  hemispheres  of  points  with  a  fixed  value 
of  B  and  L  and  (4)  some  other  geomagnetic  properties.  The 
description  includes  tables  of  the  geomagnetic  coefficients  used 
and  a  description  of  the  codes  utilized. 

The  maps  have  been  carefully  checked.  Because  of  the  magnitude 
of  this  task  and  the  wide  distribution,  notification  of  any  detectable 
errors  which  may  have  been  overlooked  would  be  appreciated. 
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DESCRIPTIVE  MATERIAL 


I  INTRODUCTION 

The  earth's  magnetic  field  has  been  studied  since  the  dawn  of  modern  physics. 
There  are  authoritative  descriptions  of  the  older  work  in  references  (1)  and  (2). 
Naturally  enought  the  earlier  students  were  interested  only  in  the  field  on  the 
surface  of  the  earth--field  which  affected  the  compass.  This  interest  led  Gauss 
to  devise  the  method  of  expansion  in  spherical  harmonics  which  has  dominated 
descriptions  of  the  geomagnetic  field  for  a  century  and  a  half. 

More  recently,  inspired  by  Birkelund's  terrella  experiments  which  so 
spectacularly  simulate  auroral  displays,  geophysicists  began  to  investigate  the 
causes  of  the  aurora  and  so  initiated  the  study  of  the  behavior  of  charged  particles 
in  space.  There  is  an  excellent  description  of  the  older  work  on  this  problem  in 
reference  (3)  and  the  newer  ideas  are  collected  in  reference  (4). 

Astrophysicists  soon  recognized  the  importance  of  this  work  to  problems  of 
stellar  evolution  and  related  astronomical  phenomena.  The  background  to  modern 
work  in  this  direction  is  well  summarized  in  reference  (5).  More  recently  still 
it  has  been  widely  recognized  that  these  effects,  originally  observed  in  an 
astronomic  context,  could  be  of  great  use  in  practical  applications  of  fusion  processes. 
The  literature  on  plasma  physics  is  already  immense  and  mostly  not  too  relevant 
to  geomagnetic  phenomena  but  some  work  which  originated  in  plasma  physics  is  of 
great  importance  (for  example,  references  (6)  and  (19). 

Finally,  the  work  of  Van  Allen  and  his  collegues  (see,  for  example,  reference 
(7)  has  clarified  the  physical  picture  involved  in  radiation  belt  and  shell  formation 
and  placed  the  mathematical  models  of  charged  particles  in  space  on  a  firm 
observational  basis. 


II  PHYSICAL  BACKGROUND 


For  an  adequate  description  of  the  graphs  in  this  report  there  are 
three  separate  stages  in  the  development  of  the  subject  which  must  be 
considered  -  first,  the  basic  geometry  of  the  earth's  geomagnetic  field;  second, 
the  motion  of  charged  particles  in  that  field;  and  third,  those  special  properties 
of  the  geomagnetic  field  of  particular  interest  to  studies  of  particle  motion. 
Each  of  these  is  discussed  separately  below. 


The  Geomagnetic  Field 


The  geomagnetic  field  unquestionably  exists  but  its  fundamental  cause 
has  not  been  satisfactorily  explained.  Paleomagnetic  studies  (9)  indicate 
that  it  has  been  present  since  the  beginning  of  geologic  history.  In  die 
absense  of  any  fundamental  theoretical  derivation  it  has  been  necessary  to 
determine  the  geomagnetic  field  strictly  from  observational  data.  This  is 
permissible  because  good  observations  exist  over  a  wide  range  of  time  and 
location. 

The  geomagnetic  induction  B  is  a  vector  quantity  which  is  conventionally 

-Hi  — ♦ 

measured  in  a  variety  of  ways.  B  is  related  to  the  force  F  exerted  on  a  particle 
of  charge  e  moving  with  velocity  u  by 


F  =  eu®  B  . 

The  geomagnetic  field  H  is  equal  to  B  divided  by  the  permeability 
of  free  space  which  is  equal  to  one  in  electromagnetic  units.  Thus 
in  electromagnetic  units  H  and  B  are  numerically  equal  but  physically  distinct 
fields.  In  this  discussion  the  three  magnetic  induction  components 


are  used  where  r,  0,  and  pare  spherical  coordinates  with  respect  to  the  earth's 
center.  The  polar  axis  will  be  identified  with  the  geographic  axis  in  such 
a  way  that  0=0  at  the  north  pole  (thus  0  measures  south  co-latitude  rather  than 


VI 


*  * 


latitude).  The  azimuthal  coordinate  tfi  is  taken  as  geographic  longitude,  east 
of  Greenwich  meridian.  The  earth's  radius  is  symbolized  by  R  and  was  taken 
to  be  6371.  2  kilometers  in  these  calculations. 

Charts  of  the  magnetic  field  at  different  epochs  are  available  from 
various  agencies.  The  United  States  Naval  Hydrographic  Office  has  under¬ 
taken  to  produce  up-to-date  compilations  at  regular  intervals.  The  data  on 
these  charts  are  used  as  input  to  analyses  described  below.  Other  agencies 
(especially  in  the  Soviet  Union)  have  also  produced  up-to-date  charts  of  the 


NORTH 


Spherical  Coordinate  System 


magnetic  field  and  theee  chart*  have  sometime*  been  used  to  modify  the 
standard  navy  charts.  Satellite  information  has  also  been  used  to  modify 
surface  observations  (see  reference  (10)). 

The  classical  Gaussian  analysis  of  the  geomagnetic  field  requires  the  use 
of  some  rather  ponderous  classical  mathematical  analysis  since  it  works  with 
expansions  in  series  of  spherical  surface  harmonics.  These  classical  mathemat¬ 
ical  techniques  are  presented  in  various  manners  in  many  different  sources. 

The  discussion  in  reference  (11)  is  rigorous  and  complete  and  provides  an 
excellent  introduction  to  the  literature.  A  more  conventional  description 
is  given  in  reference  (12). 


The  spherical  surface  harmonics  are  a  double  series  of  functions  defined 

on  the  surface  of  a  unit  sphere  with  coordinates  6  and  <p 

Gttiti  *  cos(mtp)P™(cos0)  , 

H  =  sin(mtp)  Pm(cos  0)  , 
mn  n 

where  P™  are  associated  Legendre  polynomials  as  described  below;  n  ranges 
from  zero  to  infinity  and  m  from  0  to  n.  It  should  be  noted  that  H  =  0  for 
all  n.  These  functions  are  all  mutually  orthogonal  over  the  surface  of  the 
unit  sphere,  that  is  for  all  m, m',n  and  n> 

faesine  JacpG  H  =  o, 

mn  m  n 


fd0  8in0  TScpG  G  .  ,  =  0  if  m'/m  or  nVn, 
J  1  mn  m’n' 


JaOsin©  JacpH  H  .  =  0  if  m'^m  or  n’^n, 
mn  m  n 

Thus  any  function  on  the  surface  of  the  sphere  can  be  expressed  as  the  sum  of 
series  in  terms  of  these  functions 


+  E,  h  H  ) 
mn  m=l  mn  mn 


f(0,cp)  =  2  (I  g  G 
n=o  m=o  6mn  r 

and  this  expression  is  unique  (with  the  usual  restrictions  on  discontinuities). 

As  is  shown  below  the  coefficients  g  and  h  are  easily  obtained  by 

mn  mn 

integration. 


The  associated  Legendre  polynomiala  ere  defined  in  many  different  ways. 

I 

Perhaps  the  most  convenient  form  is  the  Rodrigues  formula 


_m ,  «,  ,  ,n  (n-m).'  .  m  „  , 
Pr  (cos  9)  =  (-)  2n,  sin  9  ( 


9  ,m+n  .  n  . 
t)  sin  0. 


9  cos 


The  normalization  factor  (n-m)!  /2n!  is  that  of  reference  (13)  which  has  become 
conventional  in  a  large  part  of  geophysical  literature.  A  factor  of  l/2nn!  is  used 
in  references  (11)  (also  with  a  (-)m  multiplier)  and  (12)  and  most  standard 
references.  A  change  in  normalization  does  not  effect  the  basic  procedure 
but  it  does  change  the  values  of  the  geomagnetic  coefficients  gmn  and  hmn; 
for  example,  the  more  usual  normalization  would  lead  to  coefficients  which  are 


2n  nj  (n-m).' 

2n.' 

times  the  conventional  geomagnetic  coefficients. 


The  coefficients  g  and  h  can  be  obtained  in  the  standard  fashion 
mn  mn 


Jaesino  Jacpf(e.cp)  g  *g  raesinefopG8 

mn  mn  mn 


so  that 


2  it  (n+m).1  2n.' a 

®mn  2n+ 1  (n-m).'3  ,an  l8 

c*  n. 


JaOsinO  J  9cpGmnf(0,  V) 

®mn  2ir  (n+m)i  2n.' a 

2n+l  (n-m).'3  2ann,a 


JaesinefacpH  f(e,cp) 

.  _  _ mn _ 

mn  2ir  (n+m).1  2n.' 8 

2n+l  (n-m).'3  2ann,a 


The  particular  application  of  these  mathematical  results  which  is  of 
interest  here  is  the  case  where  f(9,^)  is  one  of  the  geomagnetic  quantities. 


I 


B  *  -  ,  S  E  m(-g  H  +  h  G  )  , 

Cp  Riino  n=o  m*o  mn  mn  mn  mn 

•o  that 

_  .  _  JaesineJacpG  bm 

R  ting  J  ^  mn  cp 

®mn  S  m  2ir  (n+m)l  2n.' a  * 

2n+l  (n-m)J 3  2ann.» 

and  so  on.  There  is  a  similar  approach  based  on  B«.  For  B_  as  will  be  dis- 
cussed  below  it  is  necessary  to  introduce  a  radial  dependence. 

The  analysis  just  described  is  an  excellent  means  of  studying  the 
geomagnetic  field  and  its  changes  so  long  as  attention  is  confined  to  the 
surface  of  the  earth.  In  order  to  describe  the  behavior  of  charged 
particles  in  space  it  is  necessary  to  extend  this  knowledge  above  the 
earth's  surface.  For  this  purpose,  some  additional  assumptions  are  required 
because  the  analysis  given  above  is  valid  for  any  function  whatsoever  defined 
over  the  earth's  surface.  If  there  were  no  geomagnetic  potential  the  only  defect 
in  the  surface  analysis  would  be  that  the  coefficients  derived  from  Bg  would 
not  equal  those  from  B^,  Since  Maxwell's  equations  virtually  insure  the 
existence  of  a  potential  this  is  a  priori  unlikely  to  occur  and  experience 
shows  that  the  minor  differences  which  do  arise  are  well  within  the  obser¬ 
vational  uncertainties. 

In  order  to  make  an  extension  of  the  analysis  into  space  (or  to  evaluate 
Br)  it  is  necessary  to  assume  that  the  values  on  the  surface  are  boundary 
values  of  some  field  which  fills  space  and  whose  behavior  is  known.  If  the 
field  is  constant  and  there  are  no  magnetic  sources  in  space  then  the 
potential  must,  by  Maxwell's  equations,  be  magnetostatic  and  satisfy 
Laplace's  equation 


0  =  hlrahv)  +  7iare^TV  + 


The  solution*  of  Laplace's  equation  (and  related  material)  are  well  known 
as  the  subject  matter  of  potential  theory  (see,  for  example,  reference  (14)). 
One  of  the  theorems  of  potential  theory  states  that  the  potential  is  uniquely 
determined  by  its  boundary  values. 


The  magnetic  field  is  definitely  not  constant.  However,  die  hope  is 
that  it  can  be  assumed  (a)  to  be  approximately  constant  (an  adiabatic 
constant  in  the  sense  used  in  the  discussion  below  on  the  behavior  of 
charged  particles)  with  respect  to  very  long  term  changes  and  (b)  to 
smooth  over  short  term  changes.  The  smoothing  has  not  been  very  carefully 
standardized- -"it  is  usual  to  take  an  annual  mean  based  on  quiet  days  as 
the  value  of  the  main  field  at  any  station"  (reference  (2)).  Many  of  the 

irregularities  arise  because  the  second  assumption - no  sources  in  space - 

is  also  invalid.  Ionospheric  currents  are  easily  detectable  but  these  alone 
do  not  cause  all  the  variations;  some  sources  of  irregularity  have  never 
been  explained.  It  is  the  hope  that  smoothing  will  permit  the  discussion 
of  any  adiabatic  constant  geomagnetic  field.  In  the  absence  of  a  definitive 
treatment  of  this  problem  it  seems  best  to  proceed  with  the  magnetostatic 
assumption. 


The  computer  codes  used  in  making  these  computations  apply  equally  well 
to  all  the  points  on  the  earth's  surface.  For  this  reason  the  data  has  been  extended 
into  close  proximity  to  the  magnet  and  geographic  poles.  Many  geophysicists 
believe  that  computations  based  on  the  magnetostatic  assumption  are  very  in¬ 
accurate  in  the  "auroral  zone"  (say,  above  about  50’  in  latitude)  because  this 
zone  is  subject  to  interference  from  interplanetary  magnetic  disturbances  and 
the  lines  of  force  originating  here  extend  so  far  into  space  that  they  may  not  close 
or  may  be  persistently  broken.  Until  this  question  is  settled  the  data  presented 
here  for  the  region  near  the  poles  should  be  used  with  caution. 


In  this  case  the  extension  from  boundary  value  to  potential  in  space  is 
formally  immediate  because 
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ia  the  unique  solution  of  Laplace'*  equation  which  vanishes  at  infinity  such  that 
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on  tiie  surface.  This  shows  how  can  be  used  to  determine  the  coefficients 
since 
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It  is  most  convenient  to  work  with  a  system  of  units  which  measure  r  in 
earth  radii  so  that  R=1  in  all  the  equations  given  above. 


It  is  easy  to  conceive  how  the  formulas  given  could  be  embodied  in 
a  computer  program.  There  is  such  a  program--a  subroutine  called 
MAGNET — which  originated  with  Air  Force  Special  Weapons  Center, 


Kirtland  A.  F.  B.  This  subroutine  exists  in  several  different  versions 
depending  on  the  values  used  for  the  geomagnetic  coefficients.  Another 
form  of  this  subroutine  utilizing  identical  mathematical  procedures  but 
written  in  FAP  for  speed  and  convenience  in  changing  geomagnetic 
coefficients  has  been  prepared  at  TEMPO.  These  versions  give 
identical  (to  five  significant  figures)  results  when  used  with  die  same 
coefficients  except  at  poles  where  the  AFSWC  version  fails  to  compute 
Bg  correctly. 

Most  of  the  computations  reported  here  used  a  set  of  48  coefficients 
obtained  by  D.C.  Jensen  and  J.C.  Cain  (reference  (15))  for  epoch  I960;  these 
were  compared  in  some  computations  to  a  set  of  512  (nominally,  actually  568 
were  used)  due  to  D.C.  Jensen  and  W.A.  Whitaker  (reference  (16))  for  epoch 
1955.  The  difference  is  relatively  small.  The  48  coefficients  of  Jensen  and 
Cain  are  shown  in  Table  A;  the  Jensen  and  Whitaker  coefficients  are  shown 
in  Table  B.  The  Jensen  and  Whitaker  coefficients  are  also  reproduced  in 
reference  (10);  the  version  in  Table  B  differs  slightly  but  the  difference  is 
probably  not  significant. 

Motion  of  Charged  Fhrticles 

The  motion  of  charged  paticles  in  a  magnetic  field  is  a  well-known 
phenomena.  In  particular,  all  the  experience  of  plasma  physics  is  available  as 
well  as  work  in  aurora  and  stellar  physics.  The  motion  has  been  studied  in 
detail  (see,  for  example,  reference  (17)).  Suppose  a  particle  has  relativistic 
mass  m  and  charge  e  and  it  moves  with  (vector)  velocity  u-  in  a  magnetic 

H  — « 

(induction)  field  B  then  its  acceleration  a  is  given  (in  rational  electromagnetic 
units)  by 

-*  a  -♦  -4  ^ 

a  =  —  u  ®  B 
m 

where  $  is  the  vector  outer  product.  If  B  is  uniform  of  strength  B  in  the 
direction  of  the  z-axis  then 

a  a  —  Bu 
x  m  y  , 


Table  A.  Geomagnetic  Coefficients  at  givei 
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= 

051253 

h 

= 
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Icienti  at  given  by  Jensen  and  Cain  (in  gauss) 


a  =  013381 

*3  8 

g  =  -.  024898 

83 

8  =  -.021795 

8  4 

g  =  -.019447 

3  6 

g  =  +.  003212 

3  6 


g  =  -.006496 

33 

g  =  +.007008 

*3  4 

g  =  -.000608 

36 
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36 
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56 
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33 
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36 
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33 
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Table  B.  Geomagnetic  Coefficients  as  given  by  Jensen  an 


GOO* 
Got* 
G  0  2  « 
G  0  3  * 
GOA* 
6  0  5  * 
G  0  6  * 
G  0  7  * 
GOB* 
G  0  9  * 
G  010  * 
G  Oil  * 
G  012  ■ 
G  013  « 
G  014  * 
G  015  * 
G  016  « 
G  017  * 
G  018  ■ 
G  019  » 
G  020  * 
G  021  * 
G  022  - 
G  023  « 
G  024  * 


-1 .8106875 7F-04 
3.03954965E-01 
2.331 168802-02 
-2. 946 138262-02 
-3.849258362-02 
1  .401941692-02 
-5.533761382-04 
-2.827916892-02 
3.967582852-02 
-5.802752202-02 
7.160703242-02 
-1.427937872-01 
6.652315972-02 
2.599230622-02 
-3.899920282-01 
4.546254102-01 
-7.676235202-01 
1.137682542+00 
6.710643892-01 
-4.416756692-01 
2.551448352+00 
9.111953022-01 
6.299493192+00 
1  .805620432+01 
-1.432658732+01 


G 

2  2*  -1.443625172-02 

H 

2  2* 

-.39770 

G 

2  3*  -2.181193712-02 

H 

2  3= 

-.40692 

G 

2  4*  -2.307357252-02 

H 

2  4* 

.87596 

G 

2  5*  -9.902658702-03 

H 

2  5= 

-.68792 

G 

2  6*  -1.336176842-02 

H 

2  6* 

-. 19722 

G 

2  7=  2.803191782-02 

H 

2  7* 

-.31370 

G 

2  8*  -5.549076142-02 

H 

2  8= 

-.50937 

G 

2  9“  5.794219852-02 

H 

2  9= 

-.23562 

G 

210*  -9.570428132-02 

H 

210* 

-.44653 

G 

211*  1.101841022-01 

H 

21  1  = 

-.26703 

G 

212*  -3.078025432-01 

H 
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-. 17468 

G 

213*  3.1 5454930E-0 1 

H 

213= 
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G 

214*  -7,381838682-01 

H 

214  = 

. 19509 

G 

215*  1.297836942+00 

H 

215= 

-.59673 

G 

216*  -1.594708932+00 

H 

216* 

•  17438 

G 

217*  6.869602862+00 

H 

217* 

-.40989 

G 

218*  -3.885785492+00 

H 

218* 

.44624 

G 

219=  -4.337150572-01 

H 

219= 

.25762 

G 

220*  —4 ,51 282966F+00 

H 

220  = 

.52052 

G 

221*  4.814281402+00 

H 

221  = 

-.29131 

G 

222=  -1.027582932+01 

H 

222= 

. 12025 

G 

223=  1 .54551627E+01 

H 

223* 

-.90117 

G 

224*  —2.61 5093472+0 1 

H 

224  = 

.28835 

G 

1  1=  1.95516829E-02 

H 
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H 
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1  2 
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G 
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G 
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1  3 
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H 

3  5* 
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G 

1  4*  —4,292285 1 4E-02 

H 

1  4 

S 

—.42734 1 04E— 02 

G 

3  6  = 

1 , 83292 1 85E-02 

H 

3  6* 

.34563 

G 

1  5*  -3.92501 289E-02 

H 

1  5 

X 

,60 30788 3E-03 

G 

3  7* 

1 . 16323982E-02 

H 

3  7* 

-.92663 

G 

1  6*  —1.51 51 9343E-02 

H 

1  6 

= 

•69540757E-02 

G 

3  8s 

-1 .21 336505E-02 

H 

3  8  = 

.13198 

G 

1  7*  9.886301 16E-04 

H 

1  7 

X 

•21372991E-02 

G 

3  9s 

6.200174  39E-02 

H 

3  9* 

.44506 

G 

1  8=  —2,6065231 6E-02 

H 

1  8 

r 

.2230874  IE-01 

G 

310s 

-1 .04773767E-02 

H 

310* 

-. 17950 

6 

1  9*  -1 .3671 1310E-02 

H 

1  9 

s 

— .38442398E-0 1 

G 

311s 

1 .69627547E— 01 

H 

31  1* 

. 1 3344 

G 

HO*  7.58558434E-02 

H 

1  10 

X 

.47657927E-0 1 

G 

312  = 

— 3 • 3358B099E-02 

H 

312* 

-. 13154 

G 

HI*  —2 . 36752740E-0 1 

H 

1  1  1 

z 

-. 13708930E+00 

G 

313  = 

3.93287235E-0 1 

H 

313  = 

.27327 

G 

112*  3.28001 598E-0 1 

H 

112 

z 

, 87569359E-0 1 

G 

314  = 

8.75455499E-02 

H 

314* 
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G 

113=  — 4.07894486E-0 1 

H 

1  1  3 

X 

•  1 1 263226E+00 

G 

315  = 

6.29087347E-01 

H 

315* 

.72092 

G 

114*  —1.21 878576E-02 

H 

114 

X 

. 35060 196E+00 

G 

316  = 

-7.85736656E-01 

H 

316* 

-.21753 

G 

115=  9.51 3738 1 6E-0 1 

H 
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X 

. 13410333E+00 

G 

317  = 

2.54874733F+00 

H 

317* 

.38167 

G 

116*  — 1 .60000 1 32E+00 

H 

1  16 

X 

•9425354 1 E+00 

G 

318  = 

-l  .068564 95E+01 

H 

318* 

•  1 48 1 9 

G 

1  17*  1  .6254371 7E+00 

H 
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,444498532+00 

G 

319  = 

7 , 36 1 68706E+00 

H 

319  = 

-.85143 

G 

118-  1  .87724 178F-01 

H 

1  18 

X 

-. 16294888E+00 

G 

320  = 

-9,39381 993E+00 

H 

320* 

•  13051 

G 

119=  —8. 2894 3944 E—0 1 

H 

1  19 

X 

— . 22336922E+00 

G 

321  * 

2, 1 73 13221 E+0 1 

H 

321  * 

-.23914 

G 

120*  1  .3035701 1F+00 

H 

120 

X 

-. 10106994F+01 

G 

322  = 

-8.0996761 3E+00 

H 

322* 

*  10810 

G 

121=  — 2.85698953F-0 1 

H 

121 

X 

•52921404E-02 

G 

323  = 

5.67008585E+0 1 

H 

323* 

— . 639 1 0 

G 

122*  -1 .25369941E+00 

H 

122 

X 

-• 1081 3708F+01 

G 

324  = 

4.85954058E+00 

H 

324* 

-.72432 

G 

1 23*  9, 1 1831 880E+00 

H 

123 

X 

• 95924 255E +00 

G 

124*  —1.21 354425E+01 

H 

124 

X 

•25457914E+01 

Coefficients  at  given  by  Jensen  and  Whitaker 


7 E-02  H  2  2* 
IE-02  H  2  3* 
5E-02  H  2  4* 
0E-03  H  2  5* 
4E-02  H  2  6* 
BE-02  H  2  7* 
4E-02  H  2  8* 
5E-02  H  2  9* 
3E-02  H  210* 
2E-01  H  21  1  = 
3E-01  H  212* 
OE-Ol  H  213* 
8E-01  H  214* 
4E+00  H  215* 
3E+00  H  216* 
6E+00  H  217* 
9E+00  H  218* 
7E-01  H  219* 
6E+00  H  220* 
OE+OO  H  221* 
3E+01  H  222* 
7E+01  H  223* 
7E+01  H  224* 


-.39770070E-02 
-•40692630E-02 
•87596 376E-02 
-.68792437E-02 
-. 19722961E-01 
-•31370295E-02 
-•5093745 IE— 02 
-. 23562953E-0 1 
—• 44653744 E—0 1 
-•26703096E-01 
-. 17468439E+00 
-• 1 1 858796E+00 

•  195096 18E+00 
-• 59673567E+00 

•  1 743809 1E+01 
-•40989000E+01 

• 44624547E+0 1 
• 25762977E+00 
. 52052 133E+01 
—.29131 02 1 E+0 1 
• 12025600E+02 
-.901 1784 2E+01 
• 28835645E+02 


G  4  4*  -2.68995824E— 03 
G  4  5*  2. 369988772— 03 
G  4  6*  2.0794081  IE-03 
G  4  7*  1 .2691 0639E-03 
G  4  8*  -6. 1  1431 164E-03 
G  4  9*  -1 .  14259993E-04 
G  410*  -1 .781 15948E-02 
G  41  1*  2.4909281  IE-02 
G  412*  -9.04762042E— 03 
G  413*  3.46934888E-02 
G  414*  7.1 924971 3E— 02 
G  415*  1 .616391 29E— 0 1 
G  416*  1 .099875 15E+00 
G  417*  -1 .69762634E-01 
G  418*  3.46227095E— 01 
G  419*  2. 12356859E+00 
G  420*  -8 .26668B59E+00 
G  421*  2.20897472E+01 
G  422*  —4 • 82690 1 44  E+0 1 
G  423*  9.254 1 2357E+0 1 
G  424=  -1 .746859 10E+02 


H  4  4*  . 14357304E-02 
H  4  5*  . 19446669E-02 
H  4  6*  . 262334 75E— 02 
H  4  7*  -.57251 951 E-02 
H  4  8*  . 69763496E-02 
H  4  9*  1 9227442E-02 
H  410*  • 14025690E-01 
H  41 1  *  . 14587651 E-02 
H  412*  . 3206085 1 E-0 1 
H  413*  .211 171 54E+00 
H  414*  -.65312709E-01 
H  415*  • 138B3826E+01 
H  416*  -.98344058E+00 
H  417*  • 42977376E+0 1 
H  418*  -.27706360E+01 
H  419*  -.24920900E+01 
H  420*  .3760514 1E+01 
H  421*  -.24900554E+02 
H  422*  . 30410382E+02 
H  423*  -. 10742B29E+03 
H  424*  . 1 1363591E+03 


IE-03 

H 

3  3* 

. 37069506E— 03 

G 

5  5* 

1 . 90808070E— 04 

H 

5  5= 

-. 39044273E— 03 

2E-03 

H 

3  4* 

. 29573649E— 02 

G 

5  6* 

—3.4 1266379E— 04 

H 

5  6* 

-•57703760E— 03 

IE-03 

H 

3  5= 

-•82332600E-03 

G 

5  7* 

1 .80337308E-03 

H 

5  7* 

. 84937564 E— 03 

5E-02 

H 

3  6* 

. 34563800E— 02 

G 

5  8* 

— 3.72495374E— 03 

H 

5  8- 

-.451 25045E— 02 

2E-02 

H 

3  7* 

-.92663568E-02 

G 

5  9* 

6 • 540 90 4 4 OE— 03 

H 

5  9= 

-•20082999E— 02 

5E-02 

H 

3  8  = 

. 131 98289E— 0 1 

G 

510* 

-3 • 2 263989 7E— 02 

H 

510* 

•  16581 872E-01 

9E-02 

H 

3  9* 

.4450536 1 E— 03 

G 

51  1* 

3.00731 552E— 02 

H 

511* 

. 33973354E— 0 1 

7E-02 

H 

310* 

-. 1 79500 16E-01 

G 

512* 

-8.50520551E-02 

H 

512* 

.364255 15E-0 1 

7E-01 

H 

31  1* 

•  1 3344460E+00 

G 

513* 

1 .81601094E-01 

H 

513* 

.66510495E— 0 1 

9E-02 

H 

312* 

-. 131 54834E+00 

G 

514* 

-3.09801 957E— 0 1 

H 

514* 

.20757329E+00 

5E-01 

H 

313* 

. 27327400E+00 

G 

515* 

2.51311 260E-01 

H 

515* 

-.36455 1 5 1 E+00 

9E-02 

H 

314* 

-.45138622E+00 

G 

516- 

-9.52850389E— 01 

H 

516- 

. 12680572E+00 

7E-01 

H 

315* 

. 72092323E+00 

G 

517* 

1 • 24242945E+00 

H 

517* 

-.1  764 34 84 E+0 1 

6E-01 

H 

316* 

-•21753822E+00 

G 

518* 

-6.841 75873E+00 

H 

518- 

-.44801 130E+01 

3F+00 

H 

317* 

•38167326E+00 

G 

519* 

5 . 0497 3274 E +00 

H 

519- 

•  152981 84E+01 

5E+01 

H 

318* 

. 1 48 1 9237E+02 

G 

520* 

-1 .06697671E+00 

H 

520* 

.465420 1 9E-0 1 

6E+00 

H 

319* 

— .85 14301 9E+01 

G 

521  * 

9.451 54250E+00 

H 

521- 

.  14300774E+01 

3E+00 

M 

320* 

. 13051 635E+02 

G 

522* 

-6 .90077794 E-0 1 

H 

522* 

. 16019346E+02 

1  E+0 1 

H 

321* 

-.23914570E+02 

G 

523. 

1 .28977227E+01 

H 

523* 

-.  1 3329579E+02 

3E+00 

H 

322* 

•  1 08 1 0833E+02 

G 

524* 

-3.41 95052 1 E+0 1 

H 

524- 

.  1 1 937782E+03 

5E+01 

H 

323* 

-•63910954E+02 

BE+00 

H 

324* 

—• 724323 14E+0 1 

1 

t - ” 

fA 


f  1 

r  ; 

t 

( 

6 

6  6* 

3.723826832-04 

H 

6  6* 

.31 1 107482-03 

6 

6  7  * 

2*488668972-04 

H 

6  7* 

•  119671652-03 

6 

6  8* 

-4.301936362-04 

H 

6  8- 

.135081652-03 

G 

6  9. 

2.137271322-03 

H 

6  9. 

-.596239412-03 

G 

610s 

3*920026092-04 

H 

610a 

-.397748392-02 

G 

61  1  S 

3.914703012-03 

H 

611a 

.143042002-01 

G 

612a 

-3.305000082-02 

H 

612a 

-.753663752-02 

G 

613a 

1 .147177382-01 

H 

613a 

. 164231522-01 

G 

614a 

-6.679614602-02 

H 

614a 

•  190444022+00 

G 

613* 

2.632547332-01 

H 

615a 

-.181180732+00 

G 

616- 

-4.683086792-01 

H 

616a 

.773856392+00 

G 

617a 

3.935888772-01 

H 

617- 

-.931853082+00 

G 

618* 

-9.467703252-01 

H 

618« 

. 118301002+01 

G 

619* 

-2.1 53346 152+00 

H 

619* 

•453616642+00 

G 

620* 

4.708427332+00 

M 

620 - 

.121374692+01 

G 

621- 

-3.239130292+00 

H 

621* 

.518860632+00 

G 

622  * 

1 .224188702+01 

H 

6 22« 

.341686982+01 

G 

623* 

-3.414990132+00 

H 

623* 

-.224784552+01 

G 

624  a 

3.475319032+01 

H 

624* 

.217233702+02 

G 

7  7a 

-3.380343522-04 

H 

.1  7a 

.212312152-03 

G 

7  8* 

3.123390792-04 

H 

7  8* 

.40 1 952472—03 

G 

7  9* 

-1 .640817672-04 

H 

7  9* 

-.3257301 12-03 

G 

710a 

-3.543096042-03 

H 

710* 

-.446024872-03 

G 

711* 

8.741389872-03 

H 

71  1- 

-.401313472-02 

G 

712a 

-2.319403382-02 

H 

712* 

. 119207182-01 

G 

713* 

9.790403632-03 

H 

713* 

-.216463872-02 

G 

714* 

-7.823776342-03 

H 

714* 

—  .251 50394  E— 0 1 

G 

715. 

4.237933782-02 

H 

715* 

— . 471 48820E-0 1 

G 

716* 

3.140930772-01 

H 

716* 

-.70200782E-01 

G 

717- 

2.350955252-01 

H 

717* 

— .63394758E+00 

G 

718- 

-2.273621892-01 

H 

718* 

. 77351 44 2E+00 

G 

719* 

-8.359978842-01 

H 

719* 

-•22506956E+01 

G 

720* 

2.118101842-01 

H 

720= 

.6029325 1E+00 

G 

721* 

-6.798134212+00 

H 

721. 

.240652372+01 

G 

722- 

1 .39271 1642+00 

H 

722- 

•40884249E+0 1 

G 

723. 

-2.0506361 72+01 

H 

723* 

• 16366798E+02 

G 

724a 

4.661321462+00 

H 

724. 

. 1 1242851E+02 
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G 

8  8.  —7.844476282-05 

H 

8  8* 

•  786K 

6 

8  9«  3.582044482-04 

+1 

8  9* 

.26501 

G 

810*  -3.298652052-04 

H 

810* 

•31681 

G 

811*  6.182627142-04 

H 

81  1  * 

-.24794 

G 

812a  -6.381378472-03 

H 

812* 

.274 1 S 

G 

813*  -9.678571382-03 

H 

813* 

.46904 

G 

814*  -1.773398382-02 

H 

814* 

-. 14794 

G 

815*  3.600073042-02 

H 

813. 

-.25111 

G 

816*  1 .26031697E-01 

H 

816. 

— . 15436 

G 

817.  -1 .30248885E-01 

H 

817* 

-. 1063! 

G 

818*  2.658526992-01 

H 

818* 

.15451 

G 

819*  3.02515981E+00 

H 

819* 

-.27101 

G 

820*  -1.313229752+00 

H 

820* 

-♦ 1631C 

G 

821*  1  .68301 795E+00 

M 

821  « 

-.20616 

G 

822*  —6.8131 2644E+00 

H 

822» 

-.18742 

G 

823-  1 .87334432E+01 

H 

823* 

-.28089 

G 

824*  -2.007614492+01 

H 

824- 

•5184! 

G 

9  9* 

2.273980682-05 

H 

9  9- 

•  10934 

G 

910- 

4,488619272-04 

H 

910* 

.21483 

G 

911* 

—5.43971 902E-04 

H 

91  1- 

-.51823 

G 

912* 

1  .2077391 3E-06 

H 

912. 

-.10370 

G 

913* 

—4 . 1 7231 458E— 03 

H 

913. 

-. 18675 

G 

914* 

-2.4091 1657E-02 

H 

914* 

•  12843 

G 

913* 

1  .533732762-02 

H 

915* 

•43221 

G 

916* 

-3.40123063E-03 

H 

916* 

.68572 

G 

917* 

-1.365730542-01 

H 

917- 

.73340 

G 

918* 

5.690847222-01 

H 

918* 

-.13342 

G 

919* 

-3.580620052-01 

H 

919* 

•25583 

G 

920* 

-1  .4295131 12+00 

H 

920* 

.32768 

G 

921* 

3,468928132-01 

H 

921* 

-.17114 

G 

922- 

-1 .140531442+00 

H 

922- 

•26847 

G 

923* 

I .706306282+00 

H 

923- 

-.56829 

G 

924* 

5.270835882+00 

H 

924* 

-.96297 

I  *  1 


i  :  i 
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28E-05 

H  8  8* 

•  786 18141 E-04 

G1010* 

2 • 264524586-05 

H1010* 

•45036482E— 04 

48E-04 

H  8  9* 

• 2650789 IE-03 

G101  1* 

2 • 335297086—04 

HI  01 1* 

-.368105296-04 

05E-04 

H  810* 

• 3 1 687960E-03 

G1012* 

-8.91 7415986-04 

H 1 0 1  2* 

-.29591 757E-03 

14E-04 

H  81  1* 

- .2479696  IE-02 

G 1  0  1  3* 

-1.568773636-04 

H 1 0 1 3* 

.2274571 16-02 

47E-03 

H  812* 

•274 1 9954E-02 

G101 4* 

-4.351263346-03 

H1014- 

.669201466-03 

586-03 

H  813* 

• 46908004 E-02 

G1015* 

-3.511591736-03 

H 1  0 1  5* 

-• 178796076-01 

38E-02 

H  814* 

-.14796194E-01 

G 1  0  1  6* 

3.328227676-03 

H 1  0 1 6* 

•33127994E— 01 

04E-02 

H  815* 

-•251 1 1 580E-0 1 

G1017* 

-5.77071 053E-03 

HI  0 1 7* 

-.82620066E— 01 

97E-0 1 

H  816* 

-• 15439786E+00 

G1018* 

-4.546206126-02 

HI  018* 

•419958896-01 

85E-0 1 

H  817* 

10635683E+00 

G1019* 

2.35214061E-01 

H 1  0 1 9* 

.21 141 051E-01 

99E-01 

H  818* 

•  1 5451 624 E+00 

G1020* 

-5. 1 0877770E-0 1 

HI  020* 

. 96 1 97826E— 0 1 

8 1 E+00 

H  819* 

-.27101 147E+01 

G1021 * 

-2.98460230E+00 

H 1 02 1 * 

. 26562 107E+00 

75E+00 

H  820* 

-• 1 63 1 0960E+0 1 

G1022* 

-4. 120526 13E+00 

HI  022* 

.20263140E+01 

95E+00 

H  821  * 

-•2061 S759E+0 1 

G1023. 

-7.74724239E+00 

HI  023* 

•B1069597E+00 

44E+00 

32E+01 

49E+01 

H  822* 

H  823* 

H  824* 

-. 187422056+01 
-.2808974 1E+02 
•518451 20E+0 1 

G1024* 

-9.34438992E+00 

HI  024* 

•43509346E+01 

58E-05 

H 

9  9* 

• 10934553E-03 

G1 1 1  1* 

3.1 1853099E-07 

HI  1  1 1* 

. 1 99435336-04 

27E-04 

H 

910* 

.21 48528  IE-03 

61 112* 

4.61401 33 | E— 05 

H1112* 

•65685716E— 04 

D2E-04 

H 

91  1* 

-.518236506-03 

G1 1 13- 

7.00662887E-04 

Hill  3* 

—.57514 181E— 03 

1 3E-06 

H 

912* 

-. 10370259E-04 

GI  1  14* 

-1 .49666325E— 03 

Hill  4* 

♦ 19025045E-02 

58E-03 

H 

913* 

-. 1 8675 1 02E— 02 

G! 1 15* 

6.207405696-03 

H 1  1  1 5“ 

-. 380370 15E-02 

57E-02 

H 

914* 

•  128455996-01 

61116* 

-4.7881 1705E-03 

HI  1 16* 

.87781 070E— 02 

76E-02 

H 

915* 

•43221558E-02 

G1  1  1  7* 

5.05607700E— 02 

H 1  1  1  7* 

-. 1 1860078E-01 

b3E-03 

H 

916* 

. 68572 128E-01 

G11I8. 

2.964810136-02 

HI  I 18“ 

-•75516208E— 02 

54E-01 

H 

917* 

. 73340895E— 0 1 

G1  1  19. 

3.37778446E— 02 

HI  1 19* 

-.277421 10E+00 

•2E-01 

H 

918* 

-• 13342141E+00 

G1 120* 

2.490354516-01 

HI  120* 

•406448316+00 

)5E-0 I 

H 

919* 

.255830376+00 

G!  121. 

5.81 554365E-01 

HI  121- 

-.626137866+00 

1  1F+00 

H 

920- 

•  327689 12E+01 

G 1 122. 

3.698043536-01 

H 1 122* 

. 965744 15E+00 

136-01 

H 

921* 

-.171148166+01 

Gl  123- 

2.05693766E+00 

HI  123* 

-•49260746E+00 

146+00 

H 

922* 

.268470656+01 

Gl 124* 

2.26466373E+00 

HI  124* 

• 588994686+0 1 

‘8E+00 

H 

923- 

-.568294156+01 

186+00 

H 

924* 

-.962974516+01 

XIII 


i  (  :  r  t 

r  t  i 

'  ( 

G1212*  4.97031879E-06 

HI  21 2“ 

1 1532726E-04 

G 1 4 14s 

1 .63016976E-05 

H 1 4 1 4* 

-.3250' 

G1 2 1 4 ■  9.32817745E-05 

H1213* 

• 662082 18E-04 

61415s 

-2.590731 38E-04 

H 1  4  1  5* 

.4950 

G1 21 5*  -1 .50069757E-03 

H 1 2 1 4* 

-• 1 156056  IE-03 

G 1  4  1  6« 

3.40091 488E-04 

H1416* 

.4822 

G 1 2 1 6*  -3.29807943F-03 

H 1 2 1 5* 

*  132 1634 7E-02 

G 1 4  1  7« 

-1 .41947019E-03 

H 1 4 1  7* 

•  2834i 

G1 2 1 7*  -2.42404020F-02 

H 1 2 1 6* 

• 39858468E— 02 

G1418- 

5.83098143E-04 

H 1 4 1 8* 

.3215 

G 1  2 1 8  ■  -1 . 18340293E-02 

H 1 2 1  7* 

•  71 293963E-03 

G 1 4  1 9« 

— 2.624055 1 2E-03 

H 1 4 1 9* 

.2033 

G 1 2 19*  -2.97000375E-02 

H 1 2 1 8* 

•54699325E-01 

G 1 4  20* 

7, 4098733 7E— 03 

HI  420* 

.5794 

G 1 220 *  -8.72789907E-02 

H1219. 

-.45257668E-01 

G 1 42 1  ■ 

-5. 1494BB57E-02 

H 1 42 1 * 

.  1 546 

G1 221  *  -1 .24087340E-01 

HI  220  s 

•24349356E+00 

G1422* 

—2. 18633866E-02 

HI  422* 

-.9805 

G 1 222*  -3.67687169E-01 

H 1 22 1  * 

-.43559620E+00 

G 1423* 

-7.  10342449E-01 

HI 423* 

-.5823 

G1223*  3.53759694E-01 

H1222* 

•86030682E+00 

G 1424* 

6. 16288781E-01 

H 1424* 

-.6828 

6 1 224 *  -4.9583547 1E+00 

H 1 223s 

-• 402624 76E+00 

H 1 224  s 

. 17689044E+01 

G 1  3 1 3* 

1 . 53294280E-05 

H 1 3 1  3  = 

- . 38266698E-05 

G1515* 

1 .21820159E-05 

HI  515* 

.2617 

G1314* 

7 . 4 564 482 0E— 05 

H 1 3  1  4  s 

. 38828054E— 03 

G1516« 

-1 .57948726E-04 

H1516* 

.3083 

G 1  3 1 5* 

-1 .44238126E-04 

H 1  3  l  5s 

—.206 121 76E— 04 

G 1  5  l  7s 

8.58468127E— 06 

H 1 5 1 7* 

.2832 

G1316. 

— 5. 87090564 E— 04 

H 1 3 1 6= 

. 24693960F— 02 

G 1 5 1  8* 

-3. 730969 13E-04 

H 1 5 1 8  * 

.  1969 

G 1  3 1  7s 

1 .70755588E-05 

H 1  3 1 7s 

•21962900E-02 

G1519* 

—2 . 27779779E— 04 

H 1519* 

•  1024 

G 1  3 1 8s 

— 5. 34952754 E— 03 

H  1  3 1 8s 

.91 944685E— 02 

G1520* 

1 .471 55267E— 02 

HI  520* 

-.6471 

G1 319* 

—6.6732664 1 E-03 

H 1  3 1 9s 

• 2 1 47B447E-0 1 

G1521* 

2.42504528E— 02 

HI  521  * 

.5167 

G 1 320* 

8. 30221 927E-02 

H  1  320* 

-.8358894  IE-01 

G 1 522* 

1  • 4 24258 1 4 E— 01 

H 1 522* 

.9672 

G1 32 1* 

—3.26044 1 78E— 0 1 

HI  321  * 

-. 1 131 4857E+00 

G1523* 

2.49694762E— 03 

HI  523* 

.4383 

G 1 322* 

7,52621 979E-0 1 

H 1 322* 

.23812565F+00 

G 1524* 

4,5455331 2E— 02 

H 1 524* 

-.4858 

G1323* 

—8.81 0603 1 4E-0 1 

H 1 323s 

-.85692 1 65 E +00 

G 1 324* 

S.95327234E-01 

H 1 324  s 

.43386539E-H30 

1 


l  :  l 


it  I  r  i 


76E-05 

HI  4  1  4* 

38E-04 

H 1  4 1  5* 

88E-04 

H 1  4  1 6* 

19E-03 

H 1 4 1  7* 

43E-04 

H 1 4  1 8“ 

12E-03 

H 1 4  1 9* 

37E-03 

H 1 420* 

57E-02 

H 1 421* 

66E-02 

HI  422* 

49F-01 

H 1423* 

81E-01 

HI  424* 

•32505932E- 05 
•49502956E-04 
•48223895E-03 
•  28346 1 6  IE-02 
•32151694E-02 
•20337665E-01 
.5794461  IE-01 
. 15466959E-01 
.980595272-01 
» 582 3203 7E+00 
.68287732E+00 


G 

1616* 

-5. 8450337  IE-06 

H 1 6 1 6* 

• 890854 26E-05 

G 

1617* 

-1 .27970049E-04 

H 1 6 1 7* 

-.  96542 149E-04 

G 

1618* 

-8.93225408E-05 

HI  6 18* 

. 70450073E-03 

G 

1619* 

-6.86292642E-04 

H 1 6 1  9* 

• 17460934E-02 

G 

1620* 

9 . 36993876E-04 

H1620* 

•  1 956297  IE-02 

G 

1621  « 

1 .42977914E-02 

HI 62 1  * 

•23798996E-0 1 

G 

1622* 

4.30473062E-02 

H1622* 

•2337769SE-0 1 

G 

1623* 

-4.1 15I7870E-02 

HI  623* 

• 15746310E+00 

fl 

1624« 

8.54049 1  BSE— 02 

H1624* 

—.1 4426384E+00 

59E-05 

HI  5 1 5* 

26E-04 

H1516* 

27E-06 

HI  5 1 7* 

13E-04 

H1518* 

79E-04 

H 1 5 1 9* 

67E-02 

H1520* 

28E-02 

HI  521* 

14E-0I 

H 1 522* 

62F-03 

HI  523* 

I2E-02 

H 1 524* 

•261 72433E— 04 
.30831 242E-04 
•28322050E-03 
. 196955 19E-02 
• 10243259E-01 
— .64714580E— 02 
.51 676833E— 0 1 
•96725937E— 0 1 
.43839802E-01 
-.48587 1 78E+00 


G 1  7 1  7*  3.58879706E— 06 
G1718«  -8.84852993E-06 
G 1 7 1 9*  -6.748881 04E-04 
G1720*  -2.30968B66E-05 
G 1 72 1 *  1 . 10360624E-03 
G1722*  1 .92448065E-02 
G 1 723*  8.300558572-02 
G 1 724  *  -4. 15158582E— 03 


H 1 7 1 7* 

. 1 25 141 04  E— 04 

Hi  7 1 8* 

. 1 5768269E— 03 

H 1 7 1 9* 

•231 50870E-03 

HI  720* 

. 1 8252449E-02 

HI  721  = 

•44591573E-02 

H 1 722* 

-.B0043459E-02 

H1723* 

.41447697E-01 

HI  724* 

-.60527826E-02 
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and  the  particle  movei  in  die  trajectory  (assuming  it  ia  initially  at  the  origin 
with  velocity  w) 

x  »  (w  +  w  ain  Xt  -  w  coa  Xt)/X  , 
y  x  y 

y  ■  (*wx  +  w^ainXt  -  w^coa  Xt)/ X  , 
z  *  w  t , 

Z 

e 

where  the  angular  frequency  X  =  —  B.  The  projection  on  the  x-y  plane  is  a  circle 

with  its  center  at  x  =  w  /X,  y  =  -w  /X  and  radius  w,  /X  where 
y  '  x  1 


r8  +  w8 
x  v 


The  projected  motion  in  the  x-y  plane  is  therefore  a  circular  path  with  period 
2ir/X.  If  the  angular  frequency  X  is  very  large  the  particle  can  then  be  smoothed 
out  over  its  circular  path  and  viewed  as  a  ring  current  of  (electromagnetic) 
intensity  eX/2ir  and  radius  w^/X  moving  along  the  z-axis  with  constant  velocity 
w  .  For  the  geomagnetic  field  near  the  earth  B  is  on  the  order  of  0.  1  gauss 
(equals  10"6  kg/sec-coulombs)  and  since  e/m  for  electrons  is  approximately 
2x  1011  coulombs/kg  it  follows  that  X  is  about  2  x  10®  revolutions  per  sec.  Since 
X  depends  inversely  on  particle  mass  the  period  of  a  proton  is  about  2000  times 
longer  than  that  of  an  electron  but  this  is  still  only  a  few  milliseconds  so  that 
the  smoothing  into  a  ring  current  seems  to  be  a  valid  physical  approximation 
for  both  particles.  The  size  of  ring  is  quite  modest--even  if  w^  is  taken  as  the 
speed  of  light  the  radius  for  an  electron  is  only  a  few  hundred  meters.  For 
protons  the  radius  is  considerable  if  the  energy  is  large  and  the  theory  given 
here  must  break  down  for  very  high  energy  protons,  say  above  1  Bev  where 
the  radius  is  about  10  km. 

Thus  the  charged  particle  will  behave  like  a  small  current  whirl 
with  the  usual  dipole  moment  (in  inductive  units) 


M  =  fcmwf/B  , 
i  0 

directed,  of  course,  in  such  a  way  as  to  oppose  the  field. 

The  physical  picture  just  described  is  often  called  the  "guiding  center 
approximation".  The  particle  is  replaced  by  a  moving  dipole  whose  position 
describes  the  trajectory  of  the  guiding  center.  In  the  uniform  field  for 
which  the  picture  was  developed  it  is  exact  but  it  is  also  applied  to 
situations  in  which  the  field  is  not  uniform.  If  non -uniformity  is  introduced 
slowly  enough  the  motion  should  be  adiabatic  (see  reference  (18)  for  a  more 
complete  discussion;  a  really  thorough  rigorous  discussion  of  this  point  is 
still  needed),  and  the  dipole  moment  should  behave  like  a  physical  invariant. 
Since  the  field  does  no  work  on  the  particle  the  total  velocity  is  a  true 
invariant.  If  w  is  the  total  velocity  then,  under  the  adiabatic  assumption. 


(w„  is  the  velocity  along  the  lines  of  force  of  the  field  and  is  the  velocity 
"perpendicular"  to  these  lines)  so  that 


=  < 


BM 

£m 


Hence  if  B  becomes  as  large  as  "x  mw'/M  it  is  necessary  for  wl(  to  go  to  zero. 


If  a  charged  particle  is  introduced  into  space  in  the  geomagnetic  region, 
say  with  initial  velocity  w(|  parallel  to  magnetic  field  and  w^  perpendicular 
to  it,  then  it  acquires  a  magnet  moment 


M  =  4  mwj  /  B 

if  B6  is  the  initial  field  strength.  It  will  move  along  the  field  lines  as 
determined  by  the  initial  component  w^ .  In  general  the  geomagnetic  field 
is  such  that  higher  altitudes  correspond  to  smaller  fields  but  each  line  of 
force  rises  to  a  maximum  and  then  returns  to  the  surface  of  die  earth.  This 
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where 


Ba  =  Ba  +  B*  =  (4coss0  +  sin36>  (m/4irrs)a  , 
r  b 

so  that  B  is  total  field  strength.  If  the  equation  for  a  line  of  force  is  utilised 
the  field  strength  at  colatitude  9  along  a  line  of  force  is  seen  to  be 

„  (3  cos3  8  + 

B  =  "_4wLs  sin8  0 -  * 

The  derivative  of  this  expression  with  respect  to  <p  is  proportional  to 
cos  9(3  +  2  cos3  0) 

so  that  there  is  a  unique  minimum  B  at  q  -  ir/2  where  B  ='^/ 4irl/.  In  this 

situation  the  pairs  of  points  with  equal  latitudes  (with  colatitudes  0  and 
ir/2-9)  are  mirror  points.  The  equator  can  be  characterized  as  those  points 
which  mirror  to  themselves. 

There  is  no  a  priori  reason  why  there  should  be  only  a  minimum  of 
B  along  any  line  of  force  but  observation  and  calculations  have  not  uncovered 
any  locations  where  multiple  minima  occur.  If  such  multiple  minima  were 
found  it  would  mean  that  some  particles  could  go  from  mirror  point  to 
mirror  point  about  one  minimum  without  even  coming  into  contact  with 
particles  moving  around  another  minimum. 

The  oscillation  of  a  particle  between  mirror  points  is  accompanied  by 
an  additional  motion  as  the  guiding  center  drifts  across  the  magnetic  field. 
This  motion  is  rather  slow  (in  the  geomagnetic  case)  and  more  difficult  to 
study  than  the  basic  case  outlined  above.  The  basic  reasons  for  th*  drift 
arise  from  the  non-uniformity  of  the  geomagnetic  force;  other  fields  also 
contribute  to  drift  but  their  contributions  are  very  small;  a  more  detailed 
discussion  is  given  below.  Plainly  it  is  necessary  to  find  some  other 
"constant"  of  motion  to  describe  the  process.  An  adiabatic  "invariant"  I 
has  been  proposed  (by  M.  N.  Rosenbluth,  see  reference  (19)  for  this 
purpose.  It  appears  to  be  based  on  die  analogy  between  a  particle  oscillating 
between  mirror  points  and  a  pendulum -type  motion  in  a  one -dimensional 
potential  where  it  corresponds  to  the  action  integral 


J  «  $W2(E-V(x))  . 

where  E  is  energy  and  V(x)  is  potential  energy. 

Since  w3  =  ws  +  -r—  B  can  be  written  in  terms  of  energy  as 
II  ?m 

Jmw1  =  Jmw*  +  MB  , 

where  Jmw1  can  be  interpreted  as  the  total  energy  available  for  motion  and 
£mw^  as  kinetic  energy  of  motion,  it  follows  that  MB  can  be  interpreted  as 
potential  energy  available  for  motion.  Then 

j  =  $ds  A«(imws-imw3^  ) 


because  M  ■  ^mwa/Bo  if  Bq  is  the  field  at  a  mirror  point.  Since ^mw  is 
a  constant  these  can  be  dropped  and 

(where  the  integral  extends  over  arc-length  from  mirror  point  to  mirror  point) 
is  an  adiabatic  invariant  (if  J  is  such  an  invariant,  see  reference  (19)).  A 
proof  that  I  is  an  adiabatic  invariant  is  presented  in  reference  (7).  Since 
the  argument  just  presented  is  highly  heuristic  such  a  proof  is  clearly 
needed. 

The  two  quantities  B0  and  I  define  a  pair  of  lines  (each  containing  the 
mirror  points  of  the  other)  around  the  world  and  therefore  would  be  adequate 
to  describe  the  drifting  of  particles  geometrically  but  not  kinematically. 

That  is,  they  determine  where  the  particles  go  but  not  the  rate  of  movement 
For  many  purposes  die  geometry  alone  is  adequate  and  it  is  this  geometry 
which  is  presented  in  the  graphs  of  this  report.  For  die  sake  of  completeness 
a  brief  sketch  of  the  details  of  motion  follows  (based  on  Alfv^n, 
reference  (5)  see  also  reference  (20)). 


There  are  two  main  caueea  for  drift  perpendicular  to  the  magnetic  field  which 
operate  in  the  absence  of  other  field* - -a  free  dipole  with  movement  M  in  a 
non-uniform  field  will  drift  a*  it  i*  acted  upon  by  a  force 

V  (M  -  B  )  . 


In  the  present  case,  however,  the  dipole  is  not  free  since  it  is  constrained  to  move 
along  a  line  of  force  and  another  force  term  must  be  included  to  account  for  the 
centripedal  effect  of  this  constraint.  This  centrifugal  force  is  given  by: 


3  a2* 

vn  as2 


where  x  is  the  position  of  the  guiding  center.  Since  x  follows  a  line  of  force 


3x  _  B 
3s  _  B 

and,  therefore 

When  the  non-uniformity  force  and  the  constraint  force  are  combined  the  total 
effective  force  is 


.  mw?  .  s  _ 
(M  +  g  11 )  V  B 


niwf. 


B 


(B  Vfl)  B 


The  second  term  of  the  expression  acts  along  the  line  of  force  and  does  not  contribute 
to  the  drift. 

In  order  to  evaluate  the  drift  consider  a  change  of  frame  to  a  system  moving  with 
drift  velocity  v  .  In  this  new  frame  the  magnetic  field  will  appear  to  contribute 
an  electric  force  given,  to  the  first  order  approximation  of  the  relativistic  correction 
by 

E  =i viS  . 


Then  the  force  on  the  particle  has  been  changed  to 
F  +  ceE'  =  F  +  ev  §  B  ; 


and  this  can  be  made  to  vaniah  by  taking 

which  is,  therefore,  the  particle  drift  velocity. 

In  terms  of  the  force  obtained  above 

Note  especially  that  since  the  charge  e  appears  it  follows  that  positive  and  negative 
particle  will  drift  in  opposite  directions--electrons  move  eastward  and  protons  west¬ 
ward  in  the  earth1  s  field.  Detailed  calculations  of  the  drift  velocity  for  a  dipole  earth 
have  been  made  several  times  (see  reference  (20)  ).  There  it  is  shown  that  a  1  MEV 
electron  requires  about  30  minutes  to  orbit  the  earth. 


Geomagnetic  Geometry 

Since  the  geometry  of  motion  of  trapped  particles  is  governed  by  the 
quantities  B  and  I  (the  subscript  aero  at  mirror  point  field  strengths  will 
be  omitted  if  there  is  no  chance  for  confusion)  it  can  be  visualised  in  terms 
of  these  quantities.  I  is, however, somewhat  awkward  to  work  with;  for 
example,  there  does  not  seem  to  be  any  closed  form  expression  for 


DP 


=  4" 


■e0  * 

sin  9  [3  cos20+l]“  [1  - 


(3  cos3  9  +  1)^ 
(3  cos30  +  lr 


sin  0  J 


the  value  of  I  for  a  particle  mirroring  at  colatitude  0O  in  a  dipole  field. 

This  is  not  especially  objectionable  on  the  basis  of  applications  because  the 
earth  is  not  exactly  a  dipole  --  the  deviation  about  5%--but  rather  because  it 
interferes  with  intuitive  understanding.  Such  understanding  is  based, 
naturally  enough,  on  the  traditional  dipole  representation  of  the  earth  which 
has  proved  very  useful.  In  am  effort  to  preserve  the  value  of  such  under¬ 
standing  Mcllwain  (reference  (21))  proposed  another  quantity  L  be  used  instead. 
Roughly  speaking  L  is  the  radial  distance  at  the  equator  which  a  line  of  force 
would  have  if  it  were  a  dipole  line  rather  than  a  realistic  one. 


Suppose  ${0)  is  the  function  defined  by  the  integral  just  quoted  so  that 


I 


DP  " 


l*<9  . 


Since,  as  shown  earlier 


and,  therefore,  if  ft  is  the  functional  inverse  of  ♦ 

13B/H  =  LaB/H#(8)  =  LaB/H*(n(L3  B/H)) 

so  that  there  is  a  functional  relationship  between  I3B/H  and  if  B/H.  This 
can  be  inv  jrted  so  that,  for  some  function  f, 

L3  B/H  =  f(l3  B/H)  . 

Thus,  if  I  is  known  L  can  be  calculated.  Mcllwain  (reference  (’21))  gives  an 
approximate  formula  for  f  and  table  which  is  presented  in  Figure  A.  A  more 
exact  formula,  also  due  to  Mcllwain  (reference  (22))  gives 

\  =  log(-~.  .  i) 

I®  B 

as  a  function  of  (j  =  log  (-jj—  )  in  the  form 

9 

X.  =  E  a  un  , 

n=0  n 

where  the  an  are  as  shown  in  Table  C.  This  formula  improves  that  of 
reference  (21)  and  was  used  (with  one  small  change)  for  all  the  calculations 
reported  here.  There  is  a  very  recent  change  included  in  Table  C  which 
was  not  present  in  the  calculations--it  seems  that  the  difference  due  to  this 
change  should  be  quite  small  (less  than  1%)  and  then  only  in  the  immediate 
neighborhood  of  the  geomegnetic  equator. 
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Figure  A.  Magnetic  Shell  Function 


Table  C.  Coefficients  in  polynomial  to  Mcllwain's  f 


If  ps-22.  0 

a0  =  +. 30062102  a!  =+.333338 

if  -22.  0  <ms  -3.  0 

l0=  +.  62337691  ai  =+.  43432642  54= +.  1501724b- 10'1 

a3= +.  13714667- 10'*  a4  =+.  82711096- 10-4  aB= +.  32916354- 1  O'5 

ae=  +.81048663-  1  o’7  a,= +.  10066362-  10-8  ^e-  +.  83232531-  10~15 

a9=  -.  81537735-  K_13 

if  -3.  0  <ms  +3-  0 

ao=+.  6228644  a1=  +.  43352788  a8  = +.  14492441-  lo'1 

a3  =  +.  11784234-  10'3  a4=  +.  38379917-  10'4  a6=  -.  33408822-  10~5 

ag=  -.53977642- 10'6  a,=  -.  21997983- 10'?  ag= +.  23028767- lo"? 

ae=  +.  26047023-  10'9 

if  +3.  0  '<ms  +11.7 

Oo=+.  6222355  a1= +.  43510529  0*=  +.  12817956-  10_1 

a3=  +.21680398-  10'3  a4=  32077032-  10'3  a6= +- 79451 31 3.  10  * 

0^=  -.  12531932- 10'4  a,=  +-  99766148-  10'®  a8=  -.  3958306- 1 0 '7 

~9=  +.  63271665-  10'9 

if  +  11.7  <U*  23 

a0  =  +2.  0007187  04=  -.  18461796  as  =  -  12038224 

a3=  -.  67310339- 10'*  a4= +.  2170224- 10‘3  aB  =-.  38049276- 10'6 

Oe=  +.  28212095-  10'7 

if  23 


04=  -3.  0460681 


04  =  +1- 0 


The  idea  behind  using  L  instead  of  I  is  that  L  is  geometrically  simplier 
than  1  but  carries  all  the  physical  significance  of  I.  This  happens  because 
L  is  almost  constant  along  lines  of  force--McUwain  reports  a  maximum 
percentage  deviation  of  1%.  All  the  particles  which  begin  to  move  from  any 
one  point  spread  themselves  along  the  line  of  force  through  that  point.  Then 
as  they  drift  perpendicular  to  the  field  their  mirror  points  sweep  out  a  surface 
defined  by  all  the  lines  of  constant  B  and  L  passing  through  the  initial  line 
of  force.  In  general,  lines  of  force  do  not  lie  in  this  surface  and  the 
particles  moving  between  the  mirror  points  do  not  move  on  the  surface.  The 
separation  away  from  the  surface  is  small  however  and  the  concept  of 
"magnetic  shell"  is  valid  if  it  is  allowed  to  have  a  finite,  though  small  thickness. 
This  shell  is  uniquely  characterized  by  L  up  to  the  variation  of  L  along  a 
line  of  force.  Hence  L  is  called  the  "magnetic  shell  radius"  or  "magnetic  shell 
parameter".  This  parameter  is  displayed  in  detail  by  the  graphs  in  this 
report. 

The  first  set  of  maps  (1  through  38)  show  contours  of  constant  B  and  L  at 
fixed  altitudes.  These  graphs  allow  data  be  reduced  in  the  B  and  L 
coordinate  system  and  then  replotted  in  geographic  coordinate  and  conversely. 

The  next  set  of  maps  (39  through  48)  show  curves  defined  by  a  fixed  pair 
of  B  and  L  values. 

The  B-L  mirror  point  traces  show  an  unexpected  geometric  feature  which 
deserves  emphasis.  All  the  points  above  any  particular  point  on  the  earth's 
surface  belong  to  different  traces  depending  on  the  values  of  B  and  L  at 
these  points.  There  is  not  a  priori  reason  why  the  different  traces  should 
not  circle  the  world  essentially  independently  of  each  other.  In  matter  of  fact, 
however,  the  different  traces  all  pass  over  virtually  the  same  points  on  the 
surface.  What  is  even  more  startling  is  that  the  opposite  mirror  points 
to  these  points  also  show  the  same  consistency  and  "stack"  over  on  another. 

Thus,  the  geometry  can  be  visualized  by  a  set  of  irregular  cones  originating 
at  the  center  of  the  earth.  These  cones  are  paired  by  containing  all  the  mirror 


XX 


points  of  one  another  and  can  be  said  to  form  a  system  of  geomagnetic 
shell  latitudes  (these  latitudes  are  illustrated  on  map  52).  The  equatorial 
cone  is  its  own  conjugate  and  is  defined  as  the  geomagnetic  equator;  it  is 
shown  on  map  51. 

III.  DESCRIPTION  OF  MAPS 

All  of  the  maps  in  this  collection  are  plotted  on  a  straight  forward  pro¬ 
jection  of  the  earth  which  permits  extension  ot  the  poles  without  unduly 
distorting  the  remainder  of  the  earth1  s  surface.  On  the  worldwide  maps  there 
is  a  superimposed  grid  of  geomagnetic  "meridians".  These  meridians  indi¬ 
cate  the  general  trend  of  lines  of  force  as  given  by  the  Jensen  and  Cain 
coefficients.  No  individual  lines  of  force  can  be  taken  as  typical  because  they 
either  do  not  reach  high  latitudes  above  the  surface  or  extend  too  far  into 
space  for  really  trustworthy  calculations.  The  approximate  meridians  shown 
were  obtained  by  following  a  "line  of  force"  constrained  to  maintain  a  con¬ 
stant  altitude  of  one  hundred  kilometers.  That  is,  it  represents  the  path  of 
the  guiding  center  of  a  charged  particle  constrained  to  this  altitude;  so  that 
it  is  not  especially  useful  in  itself  but  it  does  represent  a  convenient  approxi¬ 
mation  to  all  the  true  lines  of  force  at  different  altitudes.  These  meridians  are 
not  exact  and  should  be  supplemented  by  precise  calculations  if  accurate 
meridians  are  required. 

The  maps  are  divided  into  sections,  (A  through  G)  each  of  which  in¬ 
cluded  a  set  of  similar  displays.  The  individual  sections  are  described  in 
detail  below.  The  programs  used  are  described  very  briefly,  more  de¬ 
tailed  descriptions  will  appear  in  a  forthcoming  TEMPO  report.  All  of  the 
programs  are  available  from  the  DASA  Computer  Library  at  General 
Electric  TEMPO. 

Section  A  Contours  of  B  and  L 

Figures  1  through  16  are  maps  showing  the  contours  of  constant  values 
of  B  (the  total  magnetic  field  intensity  in  gauss)  and  L  (the  magnetic  shell 


parameter  in  earth  radii)  at  various  fixed  altitudes  above  the  surface  of  the  earth. 
These  maps  were  obtained  from  a  modified  version  of  the  computer  program  TEST 
due  to  C.  E.  Mcllwain  which  obtains  the  value  of  B  and  L  at  any  specified  location. 
The  computations  in  this  section  were  all  based  on  the  Jensen  and  Cain  set  of  48 
geomagnetic  coefficients.  The  modifications  introduced  included  a  linear  inter¬ 
polation  scheme  to  locate  contours  and  a  direct  plotter  output  subroutine  which 
provide  input  for  a  Benson-Lehner  tape-driven  off-line  plotter. 

Section  B  Contours  of  B  Near  the  South  Atlantic  Anomaly 

Figures  17  through  32  are  maps  showing  the  contours  of  constant  values 
of  B  (the  total  magnetic  field  intensity)  in  the  South  Atlantic  region  where  the 
intensity  is  minimized.  These  detailed  maps  were  prepared  because  the  intensity 
changes  very  slowly  in  this  region  and  the  large  scale  of  the  worldwide  maps 
obscures  the  structure  of  the  B  contours.  The  L  contours  are  omitted  because 
they  show  no  special  features  in  this  region.  The  computational  methods  used 
were  identical  to  these  used  in  Section  A. 


Section  C  Contours  of  B  and  L  (Jensen  and  Whitaker  Coefficients) 

Figures  33  and  34  show  the  same  material  as  section  A  (for  a  restricted 
portion  of  the  earth's  surface)  for  two  altitudes  (0  and  500  km)  with  the  difference 
that  these  maps  were  calculated  using  the  Jensen  and  Whitaker  coefficients.  Other¬ 
wise  the  methods  are  the  same  as  those  used  in  Section  A. 

Section  0  Comparison  between  Contours  (as  calculated  using  the  Jensen 
and  Cain  Coefficients  and  as  calculated  using  the  Jensen  and 
Whitaker  Coefficients) 


Figures  35  through  38  show  cross  plots  of  the  data  in  S  ection  A  and 
that  in  Section  C  for  both  altitudes  (0  and  500  km)  and  both  quantities 
(B  and  L).  These  plots  allow  the  user  to  compare  the  differences  between 
the  two  sets  of  coefficients  used.  On  a  world-wide  scale  the  difference  is 
quite  small. 
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Section  E  Constant  B-L  Traces 

Figures  39  through  44  show  the  position  and  latitude  of  the  lines 
defined  by  a  fixed  pair  of  B  and  L  for  a  variety  of  starting  conditions. 

These  lines  were  obtained  from  a  slightly  modified  version  of  the  program 

MIRROR  due  to  R.  Pennington;  MIRROR  is,  in  turn,  a  modification  of 
the  TEST  program  used  in  Section  A  so  that  all  the  calculations  are  based 
on  the  same  basic  codes.  These  calculations  used  only  the  Jensen  and 
Cain  coefficients.  As  noted  above  the  lines  originating  over  a  single 
point  (and  their  conjugate  lines)  "stack"  one  over  another  so  that  a  single 
latitude  versus  longitude  line  represents  all  the  traces  between  100  and 
1500  kilometers  to  within  1.5°  (in  the  regions  where  they  are  above  the 
surface  of  the  earth).  Superimposed  on  the  world  map  showing  the  surface 
projection  of  the  position  is  a  plot  of  the  altitude  versus  longitude. 

Section  F  Difference  between  Constant  B-L  Traces 

Figures  45  through  48  show  the  longitude  versus  latitude  dependency 
of  the  B-L  traces  for  different  altitudes  in  much  greater  detail.  This 
material  is  a  replotting  of  the  data  presented  on  a  world -wide  scale  in 
Section  E. 

Section  G  Pairs  of  Conjugate  B-L  Traces 

Figures  49  and  50  are  also  replots  of  material  in  Section  E  and 
show  the  two  conjugate  traces  simultaneously  on  the  same  map. 

Section  H  The  Geomagnetic  Equator 

Figure  51  shows  a  map  of  the  geomagnetic  equator  ot  100  kilometers 
altitude;  the  equator  at  an  altitude  of  1500  km  is  not  materially  different. 
This  map  was  obtained  by  a  special  code  which  utilized  the  subroutine 
MAGNET  with  the  Jensen  and  Cain  coefficients.  The  equator  was  obtained 
by  searching  iteratively  at  every  longitude  step  for  that  point  at  the  set 
altitude  which  was  a  minimum  of  B  along  the  line  of  force  through  itself. 
The  lines  of  force  considered  here  were  true  lines  of  force  and  since  they 
are  not  constrained  to  a  given  altitude  it  was  necessary  to  constantly 


correct  the  calculation  to  return  to  the  set  altitude.  The  equator  as  defined 
here  is  that  point  which  is  a  mirror  point  of  itself;  each  pair  of  mirror  points 
lie  one  on  each  side  of  this  equator. 

Section  I  Approximate  Geomagnetic  Shell  Latitudes 

Figure  52  shows  a  map  with  the  latitude  versus  longitude  projections 
of  a  number  of  different  traces  of  constant  B  and  L  computed  for  different 
purposes.  All  the  computations  used  the  MIRROR  code  utilized  in  Section  E 
but  there  was  no  particular  consistency  in  the  choice  of  altitudes.  Because 
of  the  miscellaneous  nature  of  the  data  these  traces  can  only  be  considered  as 
an  approximate  definition  of  the  structure  involved  in  the  geomagnetic  latitude 
system.  The  pairs  of  conjugate  cones  are  not  indicated  on  this  display  but 
can  be  inferred  approximately  from  Figures  49  and  50. 

Section  J  Conjugate  Auroral  Grids 

Figures  53  through  60  are  a  set  of  grids  showing  the  transformation 
from  (almost)  rectangular  grids  in  the  northern  hemisphere  into  grids  in 
the  southern  hemisphere.  The  points  shown  are  at  equal  altitudes  (70  and 
100  kilometers)  at  the  two  ends  of  a  common  line  of  force.  The  data  using 
the  Jensen  and  Whitaker  were  obtained  from  computations  carried  out  by 
the  Air  Force  Special  Weapons  Center  (reference  (23)).  The  data  using 
the  Jensen  and  Cain  coefficients  were  obtained  from  a  simple  line  following 
code  utilizing  the  MAGNET  subroutine.  The  altitudes  in  question  are 
especially  interesting  because  auroral  displays  occur  at  these  altitudes  and 
the  grid  shows  where  coupled  aurora  borealis  and  aurora  australis  pairs 
will  occur  simultaneously. 
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A.  CONTOURS  OF  B  AND  L 


and  Constant  Magnetic  Shell  Parameter,  L  (Earth  Radii),  At  Altitude  100  Kilometers. 
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Figure  3.  Constant  Magnetic  Field  Intensity,  B  (Gauss),  and  Constant  Magnetic  Shell  Parai 
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gure  4.  Constant  Magentic  Field  Intensity,  B  (Gauss)^  and  Constant  Magnet 


Figure  5.  Constant  Magnetic  Field  Intensity,  B  (Gauss),  and  Constant  Magnet 


B  (Gauss),  and  Constant  Magnetic  Shell  Parameter,  L  (Earth  Radii),  At  Altitude  400  Kilometers 


Figure  9.  Constant  Magnetic  Field  Intensity,  B  (Gauss),  ond  Constant  Magnetic 


,  and  Constant  Magnetic  Shell  Parameter,  L  (Earth  Radii),  At  Altitude  800  Kilometers. 
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Figure  10..  Constant  Magnetic  Field  Intensity,  B  (Gauss),  and  Constant  Magnetic  ! 


}  (Gauss),  and  Constant  Magnetic  Shell  Parameter,  L  (Earth  Radii),  At  Altitude  900  Kilometers. 
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;  (Gauss),  and  Constant  Magnetic  Shell  Parameter,  L  (Earth  Radii),  At  Altitude  1000  Kilometer?. 


Figure  12.  Constant  Magnetic  Field  Inter 


Figure  13.  Constant  Magnetic  Field  Intensity,  B  (Gauss),  and  Constant  Magnet  it 
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Figure  14.  Constant  Mognetic  Field  Intensity,  B  (Gouss),  and  Constant  Magnetic 


t,  B  (Gauss),  and  Constant  Magnetic  Shell  Parameter,  L  (Earth  Radii),  At  Altitude  1600  Kilometers. 
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Figure  15.  Constant  Magnetic  Field  Intensity,  8  (Gauss),  and  Constant  Magneti< 
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B  (Gauss),  and  Constant  Magnetic  Shell  Parameter,  L  (Earth  Radii),  At  Altitude  1000  Kilometers 

!* 
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Figure  16-  Constant  Magnetic  Field  Intensity,  B  (Gauss),  and  Constant  Magnetic 


,  B  (Gauss),  and  Constant  Magnetic  Shall  Parameter  L  (Earth  Radii),  At 

> 


2000  Kilometers 


B.  CONTOURS  OF  B  NEAR  THE  SOUTH  ATLANTIC  ANOMALY 
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insity,  B  (Gotrn),  in  the  Region  of  the  South  Atlantic  Anomaly  at  Altitude  300  km 
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Figure  21.  Constant  Magnetic  Field  Intensity,  B  (Gauss),  in  the  Region  of  the 


m  at. 


ilom«t«n. 
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Figure  22.  Constant  Magnetic  Field  Intensity,  6  (Gauss),  in  the  Region  of  th< 
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|n*ity,  B  (Gauss),  in  the  Region  of  the  South  Atlantic  Anomaly  at  Altitude  600  Kilometers. 
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Figure  25.  Constant  Magnetic  Field  Intensity,  B  (Gauss),  in  the  Region  of 


e 


^•raity,  B  (Gauss),  in  the  Ration  of  the  South  Atlantic  Anomaly  at  Altitude  800  Kilometers. 
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Figure  26.  Constant  Magnetic  Field  Intensity,  B  (Gauss),  in  the  Region  Of  th 
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|ty,  B  (Goom),  in  tho  Region  of  tho  Sooth  Atlontic  Anomaly  at  Altitodo  1000  km 
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lity,  B  (Gouts),  in  the  Region  of  the  South  Atlantic  Anomaly  at  Altitude  1200  Kilometen. 
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lity,  B  (Gauss),  in  the  Region  of  the  South  Atlantic  Anomaly  at  Altitude  1400  km 
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Figure  30,  Constant  Magnetic  Field  Intensity,  B  (Gauss),  in  the  Region  of  the  S< 
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Figure  31.  Constant  Magnetic  Field  Intensity,  B  (Gauss),  in  the  Region  of  the 
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[mity,  B  (Gouts),  in  th«  Rag  ion  of  th#  Sooth  Atlantic  Anomaly  at  Altitoda  2000  Kilomatart. 


C.  CONTOURS  OF  B  AND  L 

(Jensen  and  Whitaker  Coefficients) 


Figure  33.  Constant  Magnetic  Field  Intensity,  B  (Gauss),  and  Constant  Magnetic 


f 


,  B  (Gauss),  and  Constant  Magnetic  Snell  Parameter,  L  (Earth  Radii),  At  Altitude  0  Kilometers. 
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B  (Gauss),  and  Constant  Magnetic  Shell  Parameter,  L  (Earth  Radii),  At  Altitude  500  Kilometers. 
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D.  COMPARISON  BETWEEN  CONTOURS 

(as  Calculated  from  the  Jensen  and  Cain  Coefficients  and  as 
Calculated  by  the  Jensen  and  Whitaker  Coefficients) 


Figure  35.  Constant  Magnetic  Field  Intensity  B, (Gauss).  Using  Jensen  an< 
Field  Intensity,  B  (Gauss),  Using  Jensen  and  Whitaker's  512  Co 


Id  Intensity  B,(Gauss).  Using  Jensen  and  Cain's  48  Coefficients  and  Constant  Magnetic 
fi,  Iking  Jensen  and  Whitaker's  512  Coeficients  at  Altitude  0  Kilometers. 


Figure  36.  Constant  Magnetic  Field  Intensity  B, (Gauss),  Using  Jensen  ar 
Field  Intensity,  B  (Gauss),  Using  Jensen  and  Whitaker's  512  ( 


■Intensity  B,(Gauss),  Using  Jensen  and  Cain's  48  Coefficients  and  Constant  Magnetic 
Using  Jensen  and  Whitaker's  512  Coefficients  at  Altitude  500  Kilometers. 
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Figure  38.  Constant  Magnetic  Shell  Parameter,  L  (Earth  Radi 

and  Constant  Magnetic  Shell  Parameter,  L  (Earth  Radii),  Us 

at  Altitude  500  Kilome 


Parameter,  L  (Earth  Radii),  Using  Jensen  and  Cain's  48  Coefficients 

•r i  l  (Earth  Radii),  Using  Jensen  and  Whitaker's  512  Coefficients 
at  Altitude  500  Kilometers. 


E.  CONSTANT  B-L  TRACES 


Figure  39.  Constant  B-L  Traces  in  the  Northern  Hemisphere  at  Initial  Latitude  17. f 
Where  B  is  the  Magnetic  Field  Intensity  (Gauss)  and  L  is  the  Ma 
Note:  Above  the  Earth's  Surface  the  Latitude  Versus  Longitude  Plot  is  Common  (within  « 


irn  Hemisphere  at  Initial  Latitude  17.5°  N,  Initial  Longitude  190°  Eat  Varic 
:ield  Intensity  (Gauss)  and  L  is  the  Magnetic  Shell  Parameter  (in  Earth  Radii) 
lus  Longitude  Plot  is  Common  (within  a  maximum  deviation  t  ]  ,5°  Latitude)  I 


all  Altitude  Considered 
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Figure  40.  Constant  B-L  Trace  in  the  Northern  Hemisphere  ot  Initial  latitude  17.5  N  In 
Where  B  is  the  Magnetic  Field  Intensity  (Gauss)  and  L  is  the  Magnetic 


Figure  41 .  Constant  B-L  Traces  in  the  Southern  Hemisphere  of  the  Conjugate  to  the  Traces  at  Ir 

Where  B  is  the  Magnetic  Field  Intensity  (Gauss)  and  L  is  the  A. 
Note:  Above  the  Earth's  Surface  the  latitude  Venus  Longitude  Plot  is  Common  ( within  a  i 


Figure  42.  Constant  B-L  Trace  in  the  Southern  Hemisphere  of  the  Conjugate  to  the  Trace  at  tail 

Where  B  is  the  Magnetic  Field  Intensity  (Gauss)  end  L  is  the  W 


of  the  Conjugate  to  the  Trace  at  Initial  Latitude  17.5°  N  Initial  Longitude  190°  E  at  Initial  Altitude  1500  km 
Field  Intensity  (Gauss)  and  L  is  the  Magnetic  Shell  Parameter  (in  Earth  Radii) 


Figure  43.  Constant  B-L  Traces  in  the  Northern  Hemisphere  at  Initial  Latituc 
Where  B  is  the  Magnetic  Field  Intensity  (Gauss)  and  L  is  ths 
Note:  Above  the  Earth's  Surface  the  Latitude  Versus  Longitude  Plot  is  Common  (withii 


lorthem  Hemisphere  at  Initial  Latitude  46°  N,  Initial  Longitude  72°  E  at  Various  Initial  Altitudes 
Ic  Field  Intensity  (Gauss)  and  L  is  the  Magnetic  Shell  Parameter  (in  Earth  Radii) 

nus  Longitude  Plot  is  Common  (within  a  maximum  deviation  oft  1.5°  Latitude)  for  all  Altitudes  Considered 
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Figure  44.  Constant  B-L  Traces  in  the  Southern  Hemisphere  of  the  Conjugate  to  the  Traces  at 
Where  B  is  the  Magnetic  Field  Intensity  (Gauss)  and  L  is  the  Magr 
Note:  Above  the  Earth's  Surface  the  Latitude  Venus  Longitude  Plot  is  Common  (within  a  maxi 


phere  of  the  Conjugate  to  the  Traces  at  Initial  Latitude  46°  N  Initial  Longitude  72°  E  at  Various  Initial  Altitudes 
leld  Intensity  (Gauss)  and  L  is  the  Magnetic  Shell  Parameter  (in  Earth  Radii) 

ongitude  Plot  is  Common  (within  a  maximum  Deviation  of  t  1 .5°  Latitude)  for  all  Altitudes  Considered 


F.  DIFFERENCES  BETWEEN  CONSTANT  B-L  TRACES 
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Figure  45.  Northern  Hemiiphere  -  Difference*  in  Latitude  Venus  Longitude  Angle) 
ondThoie  Trace*  Having  the  Some  Initial  Latitude  17.5°  N,  Longil 


es  in  Latitude  Venus  Longitude  Angles  from  the  Constant  B-L  Trace  at  Initial  Altitude  of  1500  km 
Mme  Initial  Latitude  17.5°  N,  Longitude  190°  E  but  Different  Initial  Altitudes 

45 


|n  Latitude  Versus  Longitude  Angles  from  the  Constant  B-L  Conjugate  Trace  at  Initial  Altitude 
ving  the  Same  Initial  Latitude  17.5°  N,  Longitude  190°  E  but  at  Different  Initial  Altitudes 


Kern  Hemisphere  -  Differences  in  Latitude  Versus  Longitude  Angles  I 
and  Those  Traces  Having  the  Same  Initial  Latitude  46°  N,  Longitud 


Is  in  Latitude  Venus  Longitude  Angles  from  the  Constant  B-L  Trace  at  Initial  Altitude  of  1500  km 
Same  Initial  Latitude  46°  N,  Longitude  72°  E  but  Different  Initial  Altitudes 
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Figure  48.  Southern  Hemisphere  -  Differences  in  Latitude  Versus  Longitude  Angles  I 
of  1500  km  and  Those  Conjugate  Traces  Having  the  Same  Initial  Latitude  46°  N, 


is  in  latitude  Venus  Longitude  Angles  from  the  Constant  B-l  Conjugate  Trace  at  Initial  Altitude 
laving  the  Same  Initial  Latitude  46°  N,  Longitude  72°  E  but  at  Different  Initial  Altitudes 


G.  PAIRS  OF  CONJUGATE  B-L  TRACES 


Figure  49.  Constant  B-L  Trace  in  the  Northern  Hem 
of  1500  km,  initial  Latitude  17.5°  N  and  Initial  L 
Conjugate  Trace  in  the  Southern  He 


S Constant  B-l  Trace  in  the  Northern  Hemisphere  for  Initial  Altitude 
,  Initial  latitude  17.5°  N  and  Initial  Longitude  190°  E  Also  its 
Conjugate  Trace  in  the  Southern  Hemisphere 
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Figure  50.  Constant  B-L  Trace  in  the  Northern  H« 

of  1500  km.  Initial  Latitude  46°  N  and  Initial  I 

Conjugate  Trace  in  the  Southern 


H.  THE  GEOMAGNETIC  EQUATOR 


$1 .  Position  of  the  Geomagnetic  Equator  at  100km.  Altitude 


I.  APPROXIMATE  GEOMAGNETIC  SHELL  LATITUDES 


Figure  52.  Approximate  Geomagrx 


Figure  52.  Approximate  Geomagnetic  She i I  Latitudes 
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J.  CONJUGATE  AURORAL  GRIDS 
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Figure  53.  Auroral  Conjugate  Grids  for  Pacifi 
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NOTE:  Southern  Conjugate  Grids  are  Calculated  Using  Jensen  and 

Whitoker  Coefficients  and  Using  Jensen  and  Cain  Coefficients 


).  Auroral  Conjugate  Grids  for  Pacific  Area  at  70  km 
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Figure  54.  Auroral  Conjugate  Grids  for  Pacific 


SOUTH  (degrees) 


NOTE:  Southern  Conjugate  Grids  are  Calculated  Using  Jensen  and 

Whitaker  Coefficients  and  Using  Jensen  and  Cain  Coefticients 


I  Auroral  Conjugate  Grids  for  Pacific  Area  at  100  km 
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Figure  55.  Auroral  Conjugate  Grids  for  North 


.  Auroral  Conjugate  Grids  for  North  Eurasian  Area  at  70  km 
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NORTH  (  degrees) 


Figure  56.  Auroral  Conjugate  Grids  for  North 


SOUTH  (degrees) 


NOTE:  Southern  Conjugate  Grids  ore  Calculated  Using  Jensen  and 

Whitaker  Coefficients  ond  Using  Jensen  and  Cain  Coefficients 


^  Auroral  Conjugate  Grids  for  North  Eurasian  Area  at  100  km 
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Figure  57.  Auroral  Conjugate  Grids  for  W 


SOUTH  (degrees) 


SOUTH  (degrees) 


•7.  Auroral  Conjugate  Grids  for  West  Eurasian  Area  at  70  km 
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NOTE:  Southern  Conjugate  Grids  are  Calculated  Using  Jensen 

and  Whitaker  Coefficients  and  Using  Jensen  and  Cain  Coefficients 
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Figure  59.  Auroral  Conjugate  Grids  for  East  E 


SOUTH  (dearees) 


^  Auroral  Conjugate  Grids  for  East  Eurasian  Area  at  70  km 
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